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Abstract 

The exact solution for the energy spectrum of a one-dimensional Hamiltonian 
with local two-site interactions and periodic boundary conditions is determined. 
The two-site Hamiltonians commute with the symmetry algebra given by the 
Drinfeld double D{D^) of the dihedral group D3. As such the model describes 
local interactions between non-Abelian anyons, with fusion rules given by the 
tensor product decompositions of the irreducible representations of D{D^). The 
Bethe ansatz equations which characterise the exact solution are found through 
the use of functional relations satisfied by a set of mutually commuting transfer 
matrices. 

1 Introduction 

Current interest in systems of non-Abelian anyons and associated topological phases is 
primarily motivated by applications to topological quantum computing [1, 2]. There 
are indications that non-Abelian excitations occur in fractional quantum Hall systems 
with filling fraction 5/2 [3, 4]. Fractional quantum Hall systems are necessarily chiral 
in nature, breaking time-reversal symmetry due to the presence of a strong magnetic 
field. An early example of an anyonic model with chiral edge modes was given by 
Kitaev [5]. There are also expectations that topological phases may occur in time- 
reversal invariant systems [6, 7], and models have been developed in this context 
[8,9,10]. 

To help facilitate a deeper understanding of non-Abelian anyonic systems, several 
groups have undertaken studies of models of interacting non-Abelian anyons on a 
one- dimensional lattice [11, 12, 13, 14, 15, 16, 17]. Here we continue research in the 
direction of interacting systems by developing the Bethe ansatz solution for a chiral, 
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integrable chain which was introduced in [18]. The Hamiltonian of this model was 
derived using the algebraic structure given by the Drinfeld double [19] D(D 3 ) of the 
dihedral group D 3 . Consistent braiding and fusion relations for non-Abelian any- 
onic models are found in the category of irreducible representations of D(D%), since 
D(D 3 ) is an example of a quasi-triangular Hopf algebra [1, 19]. Associated with a 
(non-trivial) quasi-triangular Hopf algebra are braiding properties for anyonic degrees 
of freedom, characterised by solutions of the Yang-Baxter equation without spectral 
parameter which are realised through representations of the universal i?-matrix of 
the algebra. The anyonic fusion rules are given by decompositions of tensor product 
representations of the Hopf algebra, which are governed by the coproduct structure. 
These fusion rules provide a means to construct interacting systems by assigning 
energies to the various possible multiplet structures [13]. In [18] an integrable Hamil- 
tonian, comprised of two-site interactions which commute with the action of D(D%), 
was derived. The Hamiltonian was constructed by means of the Quantum Inverse 
Scattering Method (QISM) [20] through a solution of the Yang-Baxter equation with 
spectral parameter. However the algebraic Bethe ansatz approach to derive the exact 
solution, which often accompanies the QISM, is problematic for this model due to 
the apparent lack of a pseudo-vacuum state. 

Here we determine the Bethe ansatz solution for this model by the method of 
functional relations satisfied by a set of mutually commuting transfer matrices, fol- 
lowing techniques developed in [21, 22, 23]. While our results have close connection 
to some others in the literature, we also find differences. It can be shown [24, 25] 
that the .D (/^-invariant solution of the Yang-Baxter equation found in [18] can be 
obtained as a limiting case of the checkerboard three-state self-dual Potts model, or 
equivalently the checkerboard three-state Fateev-Zamolodchikov model [26]. However 
the one- dimensional Hamiltonian derived from this particular limiting case does not 
belong to the class of chiral Potts chain models derived in [27], this latter class arising 
from the choice of uniform rapidities. The non-chiral limit of the chiral Potts chain 
models gives the Fateev-Zamolodchikov spin models, the Bethe ansatz solution for 
which was derived in [28] using functional relations given in [22]. The energy eigen- 
values of these Hamiltonians are functions of a set of roots of a set of Bethe ansatz 
equations. What we will show below is that the energy eigenvalues of the non-Abelian 
anyonic chain we study are generically functions of two sets of roots of a single set of 
Bethe ansatz equations. This situation is reminiscent of the Bethe ansatz results for 
particular random tiling models [29, 30]. 

2 Preliminaries 

The use of Drinfeld doubles of finite group algebras for describing anyons on a lattice 
was introduced in [1]. Here we recall some relevant facts pertaining to our investiga- 
tions. 
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The dihedral group D 3 has two generators a, r satisfying: 



3 2 2 

a — e, r = e, re = a r 



where e denotes the identity. The Drinfeld double [19] of D 3 , denoted D(D 3 ), has 
basis 

{gh*\g,he D 3 }, 

where g are the group elements and g* are dual elements. This gives an algebra of 
dimension 36. Multiplication of dual elements is given by 

g*h* = 5(g,h)g* (1) 
where 5 is the Kronecker delta function. The products h*g are computed using 



h*g = g(g- l h g y. (2) 

The algebra D(D 3 ) becomes a Hopf algebra by imposing the following coproduct, 
antipode and counit respectively: 

A(gh*) = J2 ^h)* ®gk* = J2 9k* ® g{hk~ 1 )* , 
fceG fceG 

e(gh*) = 8(h,e), Vg,heD 3 . 
The universal i?-matrix is given by 

<?eD 3 

This can easily be shown to satisfy the defining relations for a quasi-triangular Hopf 
algebra: 

KA(a) = A T {a)TZ, Va G D(G), (3) 

(A®id)ft = TZ 13 TZ 23 , (4) 

(id®A)ft = 1Z 13 K 12 , (5) 
where A T is the opposite coproduct 

A T (gh*) = J2dk*® g{k- l hY = Y,9{kh~ l y ® 

fceG fceG 

The Hopf algebra D(D 3 ) has two 1-dimensional irreducible representations (ir- 
reps), four 2-dimensional irreps and two 3-dimensional irreps. Letting denote the 
matrix with 1 in the position and zeros elsewhere, the irreps are given by [18]: 
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1-dimensional irreps 



7T(i,±) = 1, T(i,±)(t) = ±1, T(i,±)(ff*) = <%> e). 

The representation 7T(i )+ ) is known as the trivial representation and is associated 
with the D(D 3 ) action on the vacuum. It can equivalently be expressed as 

7T ( i )+) (a) = e(a) Va G D(D 3 ). 

2-dimensional irreps 

Set u to be the cube root of unity 00 = e 27 ™/ 3 . Then 

*"(2,e)(<r) = ^ , 7t (2)C) (t) = ^ ^ , 7r (2) e)(^*) = 6(g,e)I 2 

and 



7T(2,i)(c) = (^q ®S)> 7T(2,i)(T)=^ ^ , 7r (2)i ) (#*) = (7)^+%, (7 ^ 

for i = 0,l, 2. 
3-dimensional irreps 




% = 0,1,2. 



Throughout we denote the module associated with an irrep 7T( aj &) by V( a ,b)- 

Tensor product representations of D(D 3 ) can be constructed through the coprod- 
uct action A. An important aspect in the definition of a Hopf algebra is that the 
coproduct is coassociative, meaning (id <S> A) A — (A <g) id)A. This allows for a 
consistently defined action of D(D 3 ) on any £-fold tensor product space. All repre- 
sentations of D(D 3 ) are completely reducible [31]. It can be shown that the £-fold 
tensor product representation of 7T( 3)+ ) decomposes as 

' \ + 1) 7T { 3, + ) © \ (3^ - 1) 7T (3 ,-), C Odd, 

<3 £ + ) = i 3 £ - 2 (7T (2 , e) © 7r (2)0 ) © 7T (2il) © 7T (2>2) ) 

©| (3 £ - 2 + 1) 7r (1)+) © | (3 £ - 2 - 1) 7r (1) _), £ even 
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where the co-efficients above denote the multiplicities. For the model we will examine 
below, the module V(3 )+ ) will be used for the local Hilbert space for the anyonic degrees 
of freedom. The global Hilbert space W is given by the £-fold tensor product 

W = V**y (6) 

However for anyonic systems, it is necessary to distinguish the physical space of states 
W C W which is governed by a superselection rule [5]. Here we take the following 
superselection rule to define W: 

ir® c +) (a)w = e(a)w Vu> G W, a G D(D 3 ), (7) 

which is equivalent to demanding that the system of D(D 3 ) anyons with local states 
represented by V(3 j+ ) are quasi particles produced from the vacuum. Mathematically, 
W is the (3 £ ~ 2 + l)/2-dimensional space of trivial modules contained in the tensor 
product module (6) when C is even. When C is odd W is zero-dimensional. 

Finally we make the following important technical observation. Let {v\, i>2, 1*3} 
denote the basis for V(3,+) with respect to the action of the defined above, i.e. 
EjV k = SjV 1 . Then the space V(i,+) C V(3,+) ® V(3,+) is spanned by the vector 

3 

s = ^v k ®v k . (8) 
k=i 

Note that s is symmetric even though the coproduct action A is not cocommutative, 
i.e. A ^ A T . 



3 The Hamiltonian and integrability 

We will study an integrable Hamiltonian, acting on W, which reads [18] 

c-i 

H = ^2 h Kk+i) + hci (9) 
k=i 

where 

h = E ® ^(3) - ^(3) ® Km) ( 10 ) 

■yeD 3 

and the elements 7 G D 3 are written as permutations of {1,2,3}. (Recall that D 3 
is isomorphic to the permutation group S3.) Since D(D 3 ) is not cocommutative, the 
inclusion of the periodic boundary interaction hci breaks the global D(D 3 ) invariance 
of H down to D 3 . Nonetheless, we will show below that the physical space of states 
W is stable under the action of H, and that the global D(D 3 ) invariance is preserved 
when we restrict the action of H to only physical states. 
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We begin by discussing some simple properties of the Hamiltonian, which is first of 
all seen to be hermitian. Defining the permutation operator P G End (V( 3)+ ) <S> V(3,+)) 
through 

P(v ®w)=W®V, V, W G V(3 )+ ) 

it can be checked that (9) commutes with the translation operator 

T = P lc Pic-i-Pi2 = exp(iV) 

where V is the momentum operator. The chiral nature of (9) can be seen through the 
space inversion of sites i >->■ C+ l — i, which induces the maps V >->■ — "P and H >->■ —if. 
That the Hamiltonian is not invariant with respect to space inversion is a signature 
of chirality. It is also apparent that the Hamiltonian is not time-reversal invariant as 
H* = —H, where throughout we will use * to denote complex conjugation. Given 
any eigenstate |\&), the time- reversed state |\&)* has energy and momentum which 
take the negative values of those for As a result we can conclude that both the 
momentum and energy spectra are symmetrically distributed around zero. 

By construction, each of the local Hamiltonians hni + \) for % = 1, C— 1 commutes 
with the action D(D 3 ) [18]. So for any a G D(D 3 ) and w G W we have 

irf 3 c +) (a)h i(i+1) w = h i{i+1) ir® c +) (a)w 
= h i(i+1) e(a)w 
= e(a)hi( i+1) w 

which shows that W is stable under the action of h^ i+ i) for i — 1, — 1. It remains 
to consider the action of hci- We can write any w G W as 

3 

k=i 

3 

= t k <g> W fc 



fe=l 



-1 



where, as a consequence of (8), {u k }\ =1 C (V(3 i+ ))' gl£ 1 and {t k }\ =l C (V(3 )+ )) 
are both isomorphic to {v k }\ =1 C V( 3j+ ) as a basis for the module associated with the 
representation 7T( 3)+ ). Using this result we have 
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Tw = ^v h (g)t k 

k=l 

which again by comparison with (8) shows that Tw G W, and necessarily T~ x w G W. 
Now we have 

h C \w = T~ x Th C \w 
= T~ l hi 2 Tw. 
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Since W is stable under the actions of h u , T and T -1 , this establishes that W is 
stable under the action of hci- Consequently W is stable under the action of H. As 
the action of D(D 3 ) is trivial on W, through the superselection rule (7), we have that 
the action of H restricted to W is both stable and .D(.D 3 )-invariant. 

Next we recall the basic results concerning the integrability of (9). We consider 
an invertible operator R(x) G End (V( 3)+ ) (g) V( 3j+ )) which satisfies the Yang-Baxter 
equation in End (V( 3j+ ) ® V( 3)+) <g) V (3i+) ): 



Ri2(x/y)R 13 (x)R 23 (y) = Rw{y)Rvz{x)R 12 {xly). 



(11) 



The argument x of R(x) is called the spectral parameter. The explicit solution related 
to (9) is [18] 

\ 



R(x) 



where 
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= x 2 
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= x, 
















d(x) 
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- X. 











(12) 



From this solution we construct the transfer matrix 

T(x) = tr (R 0C (x)R 0(c _ 1) (x)...R 02 (x)R 01 (x)) 



(13) 



which as a result of (11) forms a commuting family for different values of the spectral 
parameters: 

[T(x),T(y)] = Vx,yeC. 
Note that T(l) = T. The Hamiltonian is defined by 



H = iT~ 



(14) 



This definition give rise to the global Hamiltonian (9) as the sum of local two-site 
Hamiltonians (10) where 



h 



. d PR(x) 
dx x 



x=l 
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By construction the global Hamiltonian satisfies 

[H, T(x)\ = 0. 

As a result, T(x) can be used as a generating function for conserved operators of H. 

In order to determine the Bethe ansatz solution of the Hamiltonian, it is necessary 
to consider an expanded set of generating functions for the conserved operators. We 
introduce the operators 

1*2 



R(x) = 



{x 2 + x + l) [R'\-x 
(x 2 + x + l)(x - l) 2 



R 



-x) 



(x 2 -x + 1) 

where t 2 denotes partial matrix transposition in the second space. We also note 

1*2 



R(x) = (l-x) 2 (x 2 +x + l) 

It follows from (11) that 

Ri2(x/y)Rn(x)R 23 (y) 
Ri2{x/y)% 3 (x)% 3 (y) 
Ri 2 {x/y% 3 {x)R 23 {y) 

Ri2{x/y)R 13 {x)% 3 {y) 
Ri2(x/y)R 13 (x)R 23 (y) = 

Explicitly we have 



R \-x) 



R2z{y)Riz{x)R 12 (x/y), 

R2z{y)Riz{x)Ri2{x/y), 
R2a(y)Ria(x)Ri2(x/y), 

R2 3 (-y)Ri3(x)R 12 (x/y). 



R(x) = 



R(x) = 
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(15) 
(16) 
(17) 
(18) 
(19) 
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Next we introduce the L-operator L(x) G End (V( 2 ,i) <S> ^(3,+)) which was con- 
structed in [25, 32] and satisfies the equation 



L 12 (x/y)L 13 (x)R 23 (y) = R 23 (y)L 13 (x)L 12 (x/y) 



(20) 



in End (V( 2 ,i) <S> V(3,+) <S> V( 3j+ )). It is convenient to express the solution for the L- 
operator as 



where 



L(x) 



L\{x) 
Lj(x) 



L\{x) L\{x) 
L\{x) L\{x) 



xX, L 2 {x) = Z, 
Z~\ L 2 2 (x) = xX' 1 



and 



: 









( 1 














1' H 





UJ 





\ 


1 


J 




V 





uj- 1 ) 



X 



Using XZ = 00 1 ZX we find the inverse to be given by 

1 



L (x) 



L 2 2 (-uj- 1 x) L\(-uj- l x) 
ix 2 1 LK-ux) L 2 2 (-uj- 1 x) 



and moreover 



[L-\x)\ 



^- 2 L(-u^x). 



1 — UJ 1 x 2 
From (22) we obtain 

L 12 (x/y)L 13 (uj~ 1 x)R 23 (y) = R 23 (y)L 13 (uj- 1 x)L 12 (x/y), 
L 12 (x/y)L 13 (ujx)R 23 (y) = R 23 (y)L l3 (ux)L l2 (x/y). 

The L-operator also satisfies [25, 32] 

r 12 (x/y)L 13 (x)L 23 (y) = L 23 (y)L 13 (x)r 12 (x/y) 



(21) 

(22) 

(23) 
(24) 

(25) 



in End (V( 2 ,i) <S> V(2,i) <S> V(3,+)) where r(x) G End (V( 2 ,i) <8> V(2,i)) is a specialisation of 
the symmetric six-vertex solution: 



r(x) 



( UJX — UJ 1 x 1 














x — x 1 UJ — UJ 1 



UJ — UJ 1 X — X 1 









\ 



ujx — uj 1 x 1 J 



Now we construct the additional transfer matrices 

t(x) = tr (L oc (x)Lo {c _ 1) (x)...L 02 (x)L 01 (x)), (26) 

T(x) = ti (Roc(x)Ro(c-i){x)...Ro2(x)R i(x)), (27) 

T(x) = tT (%c(x)% { c-i)(x)..Mo2(x)R7i(x)). (28) 



As a result of equations (11,15,16,17,18,19,20,23,24,25) it follows that the set of 
transfer matrices (13,26,27,28) are mutually commuting for all values of the spectral 
parameters. This means that this set of transfer matrices can be simultaneously diag- 
onalised and their mutual eigenvectors will be independent of the spectral parameters. 
Note that T(x),T(x),T(x) are all matrices of polynomials with real co-efficients, so 
their eigenvalues will either be real polynomials or arise as complex-conjugate pair 
polynomials for x G R. To conclude this section we prove some similarly useful 
properties for t(x): 

Proposition 3.1. The transfer matrix t(x) satisfies t(—x) = (—l) c t(x), and is self- 
adjoint for iGi 

Proof Expressing the transfer matrix as 

t(x) = L%(x) ® Lg(x) ® .... ® Ljf (x) (29) 

h,...,ic 

the first part of the proposition follows from the fact that L\(x) is linear in x while 
L l j{x) is independent of x when i ^ j. That t(x) is self-adjoint for x G R follows from 

where p(l) = 2, p{2) = 1, and as a result hermitian conjugation leaves the right hand 
side of (29) invariant. □ 

4 Functional relations and the Bethe ansatz solu- 
tion 

The key observation needed in formulating a set of functional relations for the transfer 
matrices is that L{lj~ 1 ) is singular, as can be deduced from (21). It can be verified by 
direct calculations that the nullspace of Liu^ 1 ) is three-dimensional. A specialisation 
of (20) gives 

L 12 (W 1 )L 13 (x)R 2 3(u;x) = R 2 3(ojx)L 13 (x)L 12 (uj~ 1 ) 

which shows that the three-dimensional left-nullspace of L 12 (w _1 ) is invariant under 
the right action of R 2 ^{ux)Li^{x). For the tensor product of the first two vector 
spaces labelled 1 and 2, we adopt a symmetry-adapted basis consisting of the three 
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left-nullspace spanning vectors of L(uu : ) and three vectors not in the left-nullspace. 
We then find that we can write 

/ x / x ( (co^x + l)R(x) * \ 

\ {uux — 1)R(uj L x) I 

where * denotes an expression whose precise form is not needed. Through use of (30) 
we deduce that the associated transfer matrices satisfy the functional relation 

T(ux)t(x) = f(x)T(x) + g(x)T (u^x) (31) 

where we have set f(x) = (w _1 a; + 1)^ and g(x) = (ux — l) c . Making the restriction 
iGl and taking the complex conjugate of (31) leads to a second functional relation 

Tico^x) [t(x)\* = fiuj^xyTix) + g(wx)T(wx). (32) 

Starting with (23,24) and performing similar calculations leads to another four func- 
tional relations 

Tiux^iu^x) = f(u~ 1 x)T(x)+g(uj~ 1 x)T(uj~ 1 x), 
T(u- 1 x)[t(u- 1 x)Y = f(x)f(x) +g(uj- 1 x)T(ujx), 
T(ux)t{ux) = f{ux)T(x) + g{oox)T(uj- l x), 
T(u' l x) [t{ux)\* = f(ux)T(x) + g(x)T(ux). 

Let X(x), A(x), A(x), A(x) denote the eigenvalues of t(x), T(x),T(x),T(x) respec- 
tively, when acting on the global Hilbert space W (6). Keeping in mind that t(x), 
T(x), T(x), T(x) are simultaneously diagonalisable, we obtain functional relations for 
the eigenvalues which can be conveniently expressed in matrix form 

\(uj~ 1 x) —f(ou~ 1 x) —g(u~ 1 x) 
-g(x) X(x) -f(x) \ j AU'-'-.r) | - ( (J | (33) 

—f(ux) —g{ux) X(u!x) 

[X(u~ 1 x)]* —f(x) —gi^uj^x) 
-g{ux) [X(x)Y -fiu-'x) ] [ J(ux) | = | | . (31) 
-f(ux) -g(x) [\(ux)Y 

In view of Proposition 3.1 the sets of eigenvalues for t(x) and [£(#)]* are equivalent. 
To avoid confusion we have used the notations X(x) and [A(x)]* in (33,34) to make 
clear the distinction between a particular eigenvalue of t(x) compared to a particular 
eigenvalue of [i(a;)]*. 

To find the Bethe ansatz solution we make use of the facts that the matrices 
appearing in (33,34) must have zero determinant and that their nullspaces are one- 
dimensional, which can be shown to be generically true. Making the change of variable 





/ A(x) N 
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Aico^x) 




(!) 




V W / 








t A(x) \ 






)( 


A(ux) \ 


-( 


:s) 




v AVM / 
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x i — y ojx in (33), and rearranging the rows, we obtain 

Xicj-'x) -fico-'x) -g^x) \ f 1(uj- 1 x) 
-g(x) X(x) -f(x) A(ujx) 
-f(ujx) -g{ux) \(ux) J \ A(x) 

Since the nullspace is one-dimensional, we have for some rational function M(x) 

(A(uux) 
l(x) 
A^x) 

which allows us to express the nullspace spanning vector as 

A(x) 
M(x)A(ux) 
M(x)M(ux)A(u~ 1 x) 

with M(x)M(ux)M(u~ 1 x) = 1. We may now write 





where 



A(ux) M(x)A(ux) 
f(x)z 1 (x) +g(x)z 2 (x) 



. , M(ux)A(uj- 1 x) . , A(x) 

z i\ x ) = 77 — ; , Z 2{X) = 



It follows that 



A(ux) ' ZK ' M(x)A(ux)' 

zi(x)z 2 {ux) = 1, (35) 
Z\{x)z\{uix)z 1 {ui~ l x) = 1, (36) 
z 2 (x)z 2 (ujx)z 2 (uj~ 1 x) = 1. (37) 

It can verified that (35,36,37) are sufficient conditions for the matrix in (33) to have 
zero determinant. Note that Z\(x) and z 2 (x) are rational functions where the poly- 
nomials in the numerator and denominator have the same order. The most general 
form for these functions satisfying (35,36,37) is 

N N _•■ 

f = i x-iyj a± x-iyj 

for some parameters S = {i/j}^, from which an explicit expression for \(x) is found 
to be 

N N _i 

X(x) = {oo' l x + l) c TT UX ~ Wj + (ux - l) c TT - X ~ Wj . (38) 
f = i x-iyj f = \ x-iyj 

12 



Since X(x) as given by (38) must be a polynomial, setting the residues of X(x) equal 
to zero leads to the Bethe ansatz equations 



luy r l y="TT uy r Vk (39) 

when i/j ^ 0. For yj = there is no Bethe ansatz equation since the vanishing residue 
condition is automatically satisfied. 

A feature of the Bethe ansatz equations (39) is that they admit spurious solu- 
tions which do not give valid eigenvalue expressions for X(x) when substituted into 
(38). (Spurious solutions are also found for Fateev-Zamolodchikov spin models [28]). 
Consider the case N = 1 for which the Bethe ansatz equations (39) reduce to 

tooyj — 1 \ 

1 = UJ. 



icu 1 yj + 1 



The above admits solutions 



UJ 1/2C expjjni/C) + u Wexpj-jiri/C) 

In all cases except £ = 1,2, it can be verified that substituting (40) into (38) generally 
leads to an expression which is not real, in contradiction of Proposition 3.1. 
Invoking Proposition 3.1 further we calculate from (38) 



X(x) = (-l) c X(-x) 



N , . N _! 



and 

X(x) = [X(x)Y 

. r -i^r u~ x x + iy* r JL ux + iy* . , 

= ("*+i) £ n ^ +(" ^-^n^ xgm. (42) 

It is seen that (41) and (42) will be equivalent if the set of Bethe ansatz roots satisfying 
(39) is invariant under complex conjugation. This strongly suggests that the roots of 
(39) which give valid expressions for the eigenvalues X(x) through (38) will be real or 
arise in complex conjugate pairs. From numerical solutions of (39) we find that this 
is the case, with specific examples given in Table 1. For C = 1,2,3, a complete set 
of solutions for X(x) is obtained from the Bethe ansatz equations, in agreement with 
the results of computational diagonalisation. Note that in Table 1 and hereafter we 
use a two-index notation, y aj , for the roots of (39). The first index a labels the set 
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\{x) 


Bethe roots S a 


C= 1 


—x 
2x 


Si = U 

S 2 = {y2i = -l/v^} 


£ = 2 


-x 2 -l 
-x 2 + 2 
2x 2 - 1 
2a; 2 + 2 


Si = {yii = 0} 
S 2 = {y2i = -Vs/2} 
S3 = {2/31 = 0, 2/32 = -2/a/3} 

S 4 = {2/41 = 2/ 4 * 2 = - + 


£=3 


-x 3 - 2>/3siii(27r/9)a; 
-x 3 - 2 v / 3sm(87r/9)a; 
-x 3 + 2v / 3sin(47r/9)a; 
2a; 3 

2a; 3 + 3a; 

2a; 3 - 3x 


S, = {2/11 = 1.3500, 2/12 = -0.8423} 

5 2 = {2/21 = 2/22 = 0.3889 + 0.7687i} 

53 = &31 = -0.3576, 2/32 = -1.1043} 

5 4 = {2/41 = 24 = 0.2136 + 0.6230i, 2/43 = -1.3048} 

5 5 = {2/51 = 2/5*2 = 0-6704 + 1.0129*, 2/53 = -3.9134} 
S 6 = {yei = 0.7779, 2/ 62 = 0.5077, 2/ 63 = -1.4619} 



Table 1: Expressions for the transfer matrix eigenvalues A(x) and associated solutions 
of the Bethe ansatz equations for C = 1, 2, 3. For all cases the results are exact except 
the Bethe roots with C = 3 for which numerical approximations are given. Agreement 
is found between the Bethe ansatz results and computational diagonalisation. 

to which the root belongs. The second index j enumerates the elements of each set, 
which allows us to express the sets of roots as 

S a = {Vaj}f=i- 

Next we define the quantities 

Kb{%) 

where c a b is some constant and 

N a 

3=1 

It can be verified by direct substitution and use of (38,42) that the quantities (43,44,45) 
satisfy equations (33,34). Moreover we have checked for the cases C = 1,2,3 that 
(43,44,45) reproduce the transfer matrix eigenvalues which are obtained by compu- 
tational diagonalisation, with appropriate choices for the constants c ab . For the case 
C = 3 explicit values for these constants are given in Table 2, as well as the multiplici- 
ties of the eigenvalues of A a b(a;), accounting for all 27 eigenvalues. We remark that the 
degeneracies of A ab (a;), A a b(x), A a b(x) are the same for a given c a b- The occurrence of 



= c ab U^ l x) [U^)] (43) 

= Cabtaiux) MUX)}* (44) 

= c ab Ux)[Cb(x)}* (45) 



14 



a 


b 


Cab 


multiplicity of A a b(x) 


1 


1 


-3 + 2A/3sin(27r/9) 


2 


1 


2 


-2-v/3sin(87r/9) 


2 


1 


3 


-2-s/3sin(27r/9) 


2 


2 


1 


-2-\/3siii(87r/9) 


2 


2 


2 


-3 + 2V3sin(87r/9) 


2 


2 


3 


2A/3sin(47r/9) 


2 


3 


1 


—2 a/3 sinf27r/9 N ) 


2 


3 


9 


9A/3sin('47r/9"l 

y O ollll / Cf I 


9 


Q 
O 


Q 
O 


3 _|_ 9, /Q cin fA"7r /Q"\ 

o -|- za/ o oiii^vr/ y ) 


o 
L 








o 


1 
1 


4 


5 


3 


1 


4 


6 


-3 


1 


5 


4 


3 


1 


5 


5 


3 


1 


5 


6 


-3 


1 


6 


4 


-3 


1 


6 


5 


-3 


1 


6 


6 


3 


1 



Table 2: Choices of the constants c a b in (43) for C = 3 which reproduce the eigenvalues 
A a (,(x) (not shown) found through computational diagonalisation. The sets S a which 
are used to evaluate A a i,(x) are those in Table 1. The multiplicities obtained by 
computational diagonalisation are given in the final column, the sum of which confirms 
that the full spectrum is reproduced through (43). 

these degeneracies can be understood in terms of the symmetry of the problem, which 
is given by D3 since the imposition of periodic boundary conditions breaks the global 
D(D 3 ) symmetry. (Recall that for this example the physical space of states is trivial 
since C is odd.) The three-fold tensor product of a 3-dimensional irreducible rep- 
resentation of D(D 3 ) decomposes into nine copies of 3-dimensional representations. 
Each of these representations further decomposes into a sum of 1- and 2-dimensional 
representations with respect to the D 3 subalgebra of D(D 3 ). This leads in total to 
nine 1-dimensional and nine 2-dimensional irreducible representations of D 3 , which 
are in one-to-one correspondence with the results of Table 2. 

We have also confirmed that (43,44,45) hold true for £ = 1, 2, and expect that 
these formulae will give the spectrum of A (a;) for all C with real values for c a b- In 
general, assuming that the constants c a b are real ensures that the eigenvalues of 
A a b(x), A ab (x), A ab (x) will either be real or occur as complex conjugate pairs, as 
required. An unexpected feature of these results is that in many cases the transfer 
matrix eigenvalues are dependent on two set of roots of the Bethe ansatz equations 
(39). However, not all combinations of a and b give expressions in the spectrum of 
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the transfer matrices. 

As a final comment, we indicate that the constants c a } } are not determined through 
the Bethe ansatz equations (39). Fortunately, it is not necessary to know the c a b in 
order to compute the energy spectrum of (9). From the definition (14) we obtain that 
the energy eigenvalues are given by 




x=l 



(46) 



which are independent of c a b- 



5 Conclusion 

We have calculated the expression (46), which is a function of two sets of roots of the 
Bethe ansatz equations (39), as the energy eigenvalues of the Hamiltonian (9) acting 
on the Hilbert space W. We have illustrated that the equations (39) admit spurious 
solutions, and that not all combinations of roots in (46) are needed to generate the 
spectrum. A problem to address in the future is to determine the structure of roots 
which lead to legitimate expressions for the spectrum of H on W, and furthermore 
to classify the solutions which give the spectrum of H restricted to W. It would 
also be useful to characterise the appropriate root structures which correspond to the 
ground-state and elementary excitations in this latter case. 

Another direction for future work is to extend the techniques developed here to 
chains with different boundary conditions. Integrable boundary conditions for this 
model on an open chain have been determined in [33], and it is also possible to 
construct integrable closed chains with non-local (i.e. braided) boundary interactions 
between the 1st and £th sites [34]. 
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